Gap anisotropy and universal pairing scale in a spin fluctuation model for cuprates 
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We consider the evolution of d^2_y2 pairing, mediated by nearly critical spin fluctuations, with the 
coupling strength. We show that the onset temperature for pairing, T* , smoothly evolves between 
weak and strong coupling, passing through a broad maximum at intermediate coupling. At strong 
coupling, T* is of order the magnetic exchange energy J. We argue that for all couplings, pairing 
is confined to the vicinity of the Fermi surface. We also find that thermal spin fiuctuations only 
modestly reduce T* , even at criticality, but they substantially smooth the gap anisotropy. The 
latter evolves with coupling, being the largest at weak coupling. 

PACS numbers: 74.20.Mn, 74.25.Jb, 74.72.-h 



Introduction Understanding the origin of the pseu- 
dogap in high- Tc cuprate superconductors is a key prob- 
lewX. Some argue that the pseudogap originates from 
(quasi) long range order in a non-pairing channel (tvifo- 
gap scenario)^. Others argue instead that the pseudogap 
is a phase in which fermions already form singlet pairs, 
but long-range superconducting coherence is not yet es- 
tablished (one-gap scenario)'^''*. 

The theories within the one-gap scenario can be 
broadly separated into two classes - 'strong coupling' the- 
ories which consider a doped Mott insulator, and 'weak 
coupling' theories which assume a normal metallic state 
at large dopings. It is widely accepted that the cuprates 
display a crossover from Mott-like behavior in the under- 
doped regime to Fermi liquid-like behavior in the over- 
doped regime, with the superconducting dome straddling 
these two regimes. Quantitatively, strong and weak cou- 
plings are the limits of small and large values of the di- 
mensionless coupling u, which scales as U/W, where U is 
the effective Hubbard interaction, and W the bandwidth. 
For M < 1, it is natural to assume that pairing is confined 
to near the Fermi surface (FS) and can be thought of as 
mediated by a bosonic 'glue', the most natural candidate 
being collective excitations in the spin channel, enhanced 
around Q = (7r/a,7r/a) (i.e., antiferromagnetic spin fluc- 
tuations). For ii > 1, it has been argued^ that the notion 
of a bosonic glue is meaningless since the dynamics of 
the superexchange interaction, J, occurs only for energy 
scales of order U, but this picture has been challenged 
based on simulations of the Hubbard model^ . 

In this paper, we argue that the nature of the pairing 
is similar for both small and large u. We show that the 
onset temperature for pairing, which is the pseudogap 
T* in a one-gap scenario, smoothly evolves between weak 
and strong couplings, passing through a broad maximum 
at intermediate coupling, where T* is a fraction of the 
Fermi energy (Fig. [T|). At large u, T* ^ J. Still, we find 
that even in this limit, pairing is confined to the vicinity 
of the FS. The only real difference between weak and 
strong coupling is the range of the FS momenta involved 
in pairing - for u < I, pairing comes from regions around 



the hot spots (Fermi momenta connected by Q), while for 
u > I the whole FS is involved in pairing. We show that 
T* weakly depends on the magnetic correlation length ^ 
and can easily reach 300 — 500K for u ~ 0(1). 

We also discuss the special role of static thermal fluc- 
tuations, which scatter with zero energy transfer and 
therefore act as non-magnetic impurities that are pair- 
breaking for dj;2_y2 symmetry-'- . Static fluctuations are 
particularly relevant for i^ = oo, when their contributions 
to the mass renormalization and the pairing vertex di- 
verge. Earlier studies^ suggested that T* should vanish 
for ^ = OO. We found that the effect is less drastic than 
originally thought - We found that static thermal fluctu- 
ations do reduce T* somewhat, but T* still remains finite 
for ^ — cxD. 

Finally, we discuss the angular dependence of the 
dx2-y2 gap. We find that for small u, the gap is very 
anisotropic and rapidly drops in magnitude upon devia- 
tion from the hot spots. At large u, the gap near the node 
is actually larger than the cos(fca;a) — cos(/cya) form. This 
behavior can most simply be understood from the fact 
that the anisotropy of the pairing vertex $ and the mass 
renormalization Z approximately cancel in the gap pa- 
rameter A = $/Z when Z is large. More formally, static 
thermal fluctuations tend to smooth the gap anisotropy 
in order to minimize pairbreaking. 

The model We consider an approach to pairing from 
the Fermi liquid region of large dopings. We assume that 
the strongest fermion-fermion interaction is in the spin 
channel for momentum transfers near Q. The low-energy 
physics of such a Fermi liquid is captured by a spin- 
fermion model which reduces the interaction between 
low-energy fermions to the exchange of two-particle col- 
lective modes in the spin channel. 

The spin-fermion model has been described in detail in 
earlier work^. The inputs are the Fermi velocity vp, the 
spin-fermion coupling C//2, and the static boson propa- 
gator Dq{^ = 0)- In earlier work, Dq{0) was assumed to 
have an Ornstein-Zernike (0-Z) form D~^(0) = C~^+|q~ 
Qp. We will use this form of Dq{Qi), and also a related 
formi5-i(0) = r'+|q-QP+6((fe-Va)'+(9y-Va)'), 



with 6 > 0, to model inelastic neutron scattering experi- 
ments that show that the spin fluctuations decrease faster 
with deviation from Q than the 0-Z form predicts^^. 

The effective four- fcrmion vertex r((3') = {U/2)Dg gives 
rise to fermion and boson self-energies, and to an at- 
tractive pairing interaction in the dx2_y2 channel ^^. In 
Eliashberg theory-''^^, which has been justified earlieriSii^^ 
the linearized gap equation has the form 
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where w, uj' are Matsubara frequencies, the integration 
over k' is along the FS, i^k is the density of states nor- 
malized such that :/ = 1 along the nodal direction, and 
u = 3Ua/{8'KVp) is the dimensionless coupling where 
vp is the nodal velocity. The full dynamical interaction 
Dq{n) includes the boson self-energy, which in Eliashberg 
theory reduces to the Landau damping term 
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We emphasize that the boson self-energy comes from the 
same interaction that gives rise to pairing, and T contains 
the same dimensionless coupling u. 

The overall sign in the r.h.s. of Eq. ([T]) is a consequence 
of the fact that pairing is in the spin channel. Isotropic 
s— wave pairing is not possible in this case, while d^2_y2 
pairing is favored because Afe(a;) and Afc/(w) have oppo- 
site sign for k' = k + Q, where Dq^q) has a maximuniii. 

Eq. Ill has a solution at T = T* . As the temperature 
only appears explicitly in Eq. ([1]) in the Landau damping 
term fl/T = 2TmT/r, T* ex P. Using the definition of F, 
we then obtain 



T* 



-k{u) 
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where /c(u) is a function of u and ^. 

We also consider the mass renormalization Zk{uj) = 
1 -1- Sfe(a;)/(jj, which effectively measures the strength of 
the coupling along the FS. It is given by 



uT 



Zk{uj) = 1 + Y^ jdk'vk 



signuj 



Dk-k'{u:-uj') (4) 



One can easily verify that ZkiirT) = 1 + 
{u/ (an)) f dk' i^k' Dk-k' (0) is independent of the Landau 
damping. For large u, Z ^ 1 along the entire FS; for 
small u, but large ^, it is still large near the hot spots 
k = kh, where Zkf^{nT) = 1-1- u{^/a). We verified that 
for ^ ^ a, fermions relevant for pairing have Zk > I, 
i.e., pairing in the critical regime is a strong coupling 
phenomenon for all u. 

The results The results of the numerical calculations 
for T*{u), Afe(w) and Zk{uj) for the 0-Z form of D,(0) 
are presented in Figs. [T][3] (we used 142 Matsubara fre- 
quencies). Akiuj) monotonically decreases with increas- 
ing frequency (not shown). The dependences of Akico) 




FIG. 1: (Color online) The onset temperature T* for d^2_y2 
pairing vs the dimensionless coupling u for an Ornstein- 
Zernike form for the static spin propagator near Q and various 
values of the magnetic correlation length ^, with vp/a = leV. 
The insets show where the relevant fermions are located along 
the FS. For both small and large u, pairing is confined to its 
vicinity. 
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FIG. 2: (Color online) The mass renormalization Zk{TrT) and 
the pairing gap Afe(7rT) for small u = 0.2. Panels (a) and 
(b) — the results for two different ^; panel (c) - the angular 
dependence of the gap compared to the cos{kxa) — cos{kya) 
form; panel (d) - Ak{TrT) for ^ = 4.8a with and without 
the pairbreaking contribution from the static spin fiuctua- 
tions (the term with uj — uj' in Eq. ([TJ). Without the static 
contribution, ZkinT) = 1. 



and Zk{uj) for k along the FS are similar for different w, 
so we only present the results for the lowest Matsubara 
frequency uj = ttT. 

In Fig. [l] we show the dependence of T* on the dimen- 
sionless coupling u. We see that T* initially increases 
with u, passes through a broad maximum at u ^^ 1, and 
then decreases, eventually as 1/u. 

In Fig. [2] we show the results for Zk{nT) and Ak{nT) 
for small u — 0.2. We see that ZkinT) is enhanced near 
the hot spots and then drops to near its non-interacting 
value (Z = 1), which implies that only the hot regions 
are relevant for pairing. Afe(7rT) also has a maximum at 
the hot spots, leading to a significant deviation from the 
cos(fca;a) — cos(fcj,a) form. Moreover, the slope of Ak{uj) 
near the node is smaller than this form. 
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FIG. 3: (Color online) Same as Fig. (2] but for large u = 2. 
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FIG. 4: (Color online) The onset temperature T*{u) for dif- 
ferent forms of the static susceptibility Dq{0) (see text). The 
magnitude of T*{u) increases as Dq{0) drops faster with de- 
viation from Q. 



In Fig. [3] we show the results for large u — 2. We see 
that Z is enhanced along the entire FS. The deviation 
of the gap anisotropy from the cos(A:a;a) — cos(fcya) form 
is much weaker than at small u, and the slope of Afe(a;) 
near the node is larger than this form. 

In Fig. |4] we show how T* (u) changes with the devia- 
tion of the static susceptibility from the 0-Z form. We 
see that the magnitude of T*{u) increases for b > 0, i.e., 
when -05(0) drops faster from Q than the 0-Z form. 

Analytical reasoning All these results can be under- 
stood analytically. Consider first why T* remains fi- 
nite when ^ diverges, despite the fact that thermal spin 
fluctuations are pairbreaking for a d— wave gap. We re- 
call that Eliashberg theory is a set of two coupled equa- 
tions for the pairing vertex ^{uj) and self-energy T,{uj). 
They reduce to independent equations for A(w) and 
Z{ui) by exploiting the definitions <l>(a;) = A{uj)Z{ui) and 
S(a;) = uj{Z{uj) — 1). The static thermal contribution to 
Eq. Ill) (the term with uj = lo') has the form 



a\uj\ J ^ 



Afc'(w)-hAfc(cj) 



|k-k'-Q| 



(5) 



where the terms with Afe and A^/ are the contributions 
from S and $, respectively. Each diverges at ^ = cxd, and 



taken alone, would drive T* to zero. However, the sum 
of the two remains finite because by symmetry A^+q — 
— Afc, and T* does not vanish at 1^ = 00. Note that there 
is no such cancellation for Z{ijj), which diverges at a hot 
spot when ^ = 00 (Ref. 16) 

Consider next the dependence of T* on u for ^ = 
00. At small u, pairing is confined to hot regions 
{k « kh)^ and the momentum dependence of the static 
Dk~k' can be approximated by |k;i — k'p. The inte- 
gral over k' in Eq. ([T]) can then be evaluated analyt- 
ically, and Eq. ([1]) reduces to a one-dimensional inte- 
gral equation for Afe^(a'). Simple power counting then 
yields T* ~ (r/a^)M^. A numerical solution give^ii T* « 
0.68(r/a2)u2 = 0.13{vF/a)u, i.e., /oo(u < 1) « 0.13u. 

The reduction of T*{u) at large u is peculiar to 
d— wave pairing. The argument is that, as u increases, 
T* initially also increases, and at some u the dynamic 
term in Dk-k'{^ — ^') becomes comparable to a typi- 
cal |k — k' — Qp term along the FS, which determines 
the attractive d— wave component of the static interac- 
tion. A further increase of T would make the effective 
interaction less momentum dependent and hence would 
reduce the d— wave attraction. The balance is reached 
when T* /T is a constant, i.e., when foc{u) oc 1/u. Nu- 
merically, foa{u) ^ 0.056/m for the 0-Z form of Dq{0). 
This can be re-expressed as T* ~ 0A7{vF/a)'^/U ~ 0.5 J, 
where J ~ (vp/a)'^ /U is the exchange integral of the un- 
derlying Heisenberg model at half fiUingi^. We interpret 
this as evidence that pairing of incoherent fermions in 
the spin-fermion model, and the creation of singlet pairs 
upon doping a Mott-Heisenberg insulator, describe the 
same physics from different perspectives. 

At the same time, we argue that, even at strong cou- 
pling, pairing is confined to the vicinity of the FS. To 
see this, from the integral over ek which leads to Eq. ([1]), 
we can estimate a typical |k — ki?| transverse to the FS 
from vplk — ki?| ~ luZ(uj), where uj ~ ttT* is a typical 
frequency for the pairing problem. Using the quantum 
critical form^ Z{uj) « 2uyT/w, we find that for large u, 
k - k^l ^ LLitypZ{ujtyp)/vF ^ (3/8)v^7rr*/r ~ 0.1(7r/a). 
This |k — k^l is numerically much smaller than kp ~ 
0.8(7r/a), i.e., pairing involves fermions from a narrow 
shell around the FS. For smaller u, |k — kp\/kp is even 
smaller. 

Consider next the variation of the gap along the FS. 
For small u, the anisotropy of the gap is a consequence of 
the fact that the pairing problem is confined to a region 
near the hot spots with width Sk ~ kpu. We found that 
the ratio of the slope of the gap near the node to the 
gap value at the hot spots scales as w^ for small u. At 
strong coupling, the non-confinement of A/j to hot spots 
is due to a cancellation between the anisotropics of $ and 
Z when dividing to form A. Further, expanding Eq. ([5|) 
near k' = k + Q, we find that the residual, non-singular 
pairbreaking contribution from thermal spin fluctuations 
contains d^A/d^k', i.e., it tends to make the gap more 
smooth. The extent to which this affects the shape of 
the gap can be verified numerically. In the last panels in 



Figs. [2] and [3l we show the gap variation along the FS 
with and without the pairbreaking contribution to the 
gap equation. We see that Afc near the node becomes 
larger than the cos{k,j;a) — cos(kya) form when the pair- 
breaking contribution is included. 

Finally, we consider the variation of the overall scale for 
T* with b (Fig. 2]). To understand this, we recall that for 
the 0-Z form oi Dq{0), a strong reduction of T* compared 
to the asymptotic, small u form T* w 0.13{vF/a)u was 
due to the interplay between the \lu — uj'\/T term and the 
maximum |k — k'| along the FS. Once ^^^(O) becomes 
steeper with deviation from Q (& > 0), the typical k and 
k' get closer to the hot spots, the restriction on |k — k'| 
becomes less relevant, and T* increases. For 6 < 0, the 
situation is the opposite, geometrical restrictions become 
more relevant, and T* rapidly decreases. 

Comparison with earlier studies Several groups did ex- 
tensive studies of T* within the Eliashberg theory for the 
spin-fermion modeli^iiilii^. They treated F as an inde- 
pendent parameter, i.e., they did not express it in terms 
of u. Our results agree with these studies if we use the 
same vp, u and F as they did. Monthoux and Pinesii 
studied T*{u) at small u and noticed that its slope de- 
creases as the coupling increases, consistent with Fig. [1] 
Schuttler and Normar>i^ found a saturation of T* when 
measured in units of F, but for a model in which D[q^ Vi) 
had the factorizcd form Z?i(g)Z?2(^), and a strong T de- 
pendence of F is introduced phenomenologically. Within 
the spin-fermion model, though, the T dependence of F 
is weak. Our results also agree with FLEX calculations 
for the Hubbard modeli^i^Sii^i. In these calculations, F is 
obtained self-consistently, as in our theory. For U = At 
and vp/a « 2i (u w 0.25]s our T* « 0.02t is in good 
agreement with Refs. [l9li20li2ll . This value is also in 
good agreement with two-particle self-consistent calcu- 
lations2^, dynamical cluster approximation^^, and clus- 
ter DMFT—^, which also yields that the d— wave order 
parameter scales as J at large u (Ref. \23 i) . 

Comparison with the cuprates To get T* in Kelvin, 
we use vp/a '^ leV, noting that vp is the 'bare' velocity 
as obtained in band theory. For the O-Z form of the 
static Dq{0), we obtain T* - 0.02vp/a - 200 - 250/^ for 



u ~ 0(1). For 6=1, this temperature increases to over 
350K (Fig.H]), and becomes even larger for larger b. This 
shows that a spin-mediated pairing interaction is strong 
enough to account for experimental values of T* . 

We also verified that the 'hot spot' description, which 
yields a very anisotropic d— wave gap, is only valid for 
small u, which would correspond to strongly overdoped 
cuprates. But in this case, ^ is also small, and therefore 
the enhancement of Afc near the hot spots is weakened. 
For optimal and underdoped cuprates, we find a gap close 
to the cos{kxa) — cos(fcj,a) form. This is consistent with 
photoemission^S below Tc, and with the scenario that the 
Fermi arcs above Tc appear because of thermal broaden- 
ing of the spectral function''. 

Conclusions In this paper, we analyzed how (1^2 _y2 
pairing mediated by nearly critical spin fluctuations 
varies with the coupling strength. We argued that 
the onset temperature for pairing T* « {vp/a)f^{u) 
smoothly evolves between weak and strong couplings, 
passing through a shallow maximum at w ~ 1. At large 
u, T* '^ vp/u r^ J. For all u, pairing is confined to the 
vicinity of the FS. We also argued that T* only weakly 
depends on the distance to the antiferromagnetic insta- 
bility, as singular pairbreaking contributions from static 
spin fluctuations cancel out in the gap equation. The 
residual pairbreaking terms only modestly reduce T*, 
and at the same time smooth the angular dependence 
of the gap, and at strong coupling make A^ near the 
node even larger than the cos(fca;a) — cos{kya) form. As 
Mott physics is certainly present near half filling, where 
T* reaches its largest values, our T*{u) should only be 
taken as an estimate. Still, the fact that T* is in the 
experimental range is in support of a one-gap scenario in 
which the instability at T* occurs in the particle-particle 
channel, due to interactions with spin fluctuations. 
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